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A Poincare covariant Faddeev equation is presented, which enables the simultaneous prediction 
of meson and baryon observables using the leading-order in a truncation of the Dyson-Schwinger 
equations that can systematically be improved. The solution describes a nucleon's dressed-quark 
core. The evolution of the nucleon mass with current-quark mass is discussed. A nucleon-photon 
current, which can produce nucleon form factors with realistic Q 2 -evolution, is described. Axial- 
vector diquark correlations lead to a neutron Dirac form factor that is negative, with r™ u > r™ d . 
The proton electric- magnetic form factor ratio falls with increasing Q 2 . 
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In quantum field theory a nucleon appears as a pole 
in a six-point quark Green function. The residue is 
proportional to the nucleon's Faddeev amplitude, which 
is obtained from a Poincare covariant Faddeev equa- 
tion that sums all possible exchanges and interactions 
that can take place between three dressed-quarks. A 
tractable Faddeev equation for baryons was formulated 
in Ref. [1]. Depicted in Fig.[TJ it is founded on the obser- 
vation that an interaction which describes colour-singlet 
mesons also generates quark-quark (diquark) correlations 
in the colour-3 (antitriplet) channel [2[. While diquarks 
do not appear in the strong interaction spectrum; e.g., 
Refs. the attraction between quarks in this chan- 

nel justifies a picture of baryons in which two quarks are 
always correlated as a colour-3 diquark pseudoparticle, 
and binding is effected by the iterated exchange of roles 
between the bystander and diquark-participant quarks. 

The study of baryons using the Poincare-covariant 
dressed-quark Faddeev equation sits squarely within the 
ambit of the application of Dyson-Schwinger equations 
(DSEs) in QCD. Notably, the DSEs admit a nonpertur- 
bative symmetry-preserving truncation scheme 0, 0, @ , 
which has enabled the proof of numerous exact results 
for pseudoscalar mesons, and also provides for the for- 
mulation of reliable models that can be used to illustrate 
those results and make predictions for a wide range of 
meson observables with quantifiable errors; e.g., Ref. [7[. 





FIG. 1: Covariant Faddeev equation, ip: Faddeev amplitude 
for nucleon of total momentum P = p q + pd- It expresses 
the relative momentum correlation between the dressed-quark 
and -diquarks within the nucleon. Shaded region, equation's 
kernel: single line - dressed-quark propagator, S(p), see 
Eq. dl]); T - diquark correlation amplitude, see Eq. J5]); and 
double line - diquark propagator, see Fig. [3] 



Progress was recently made toward placing the com- 
putation of nucleon properties on the same level as that 
of mesons [8|- In enabling the direct correlation of 
meson and nucleon properties via a single interaction 
kernel that preserves QCD's one-loop renormalisation 
group behaviour, it significantly complements existing 
phenomenological studies; e.g., Refs. 0, Qij. Neverthe- 
less, Ref. [8] can be refined. Herein we report two mate- 
rial steps in that direction. 

The first is to base the Faddeev equation on 
the amended one-parameter renormalisation-group- 
improved rainbow-ladder interaction introduced in 
Ref. [U|; i.e., with i = k — q, 
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with Z2 the fermion wave function renormalisation con- 
stant, F(s) = [1 - exp(-s/A 2 )]/s, A t = 1.0 GeV, r = 
e 2 - 1, 7m = 12/25, and A QC d = 0.234 GeV. In Eq. © 

C\(uj — uj(m)) ^ ^ 
m) = Co + 1+C2£ + C 2^ x = W-o, (3) 

d(t) = 0.86(1 - 0.15* + (1.50t) 2 + (2.95t) 3 ) , 

with uj(m) = 0.38 + 0.17/(1 + m/mo), where m is 
the renormalisation-group-invariant current-quark mass, 
m = 0.12 GeV, and C = 0.11, C 2 = 0.885, C 3 = 0.474. 
NB. By including dependence on m, limited aspects of 
vertex corrections are represented in the kernel. 

As explained in Ref. [ll| , the interaction in Eq. ([2|) 
is characterised by a single constant; viz., C(o;,m). It 
was determined subject to an understanding that cor- 
rections to the rainbow-ladder truncation which impact 
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upon hadron phenomena, vanish with increasing current- 
quark mass; and in connection with light-quark systems, 
and those of the physical qualities of the pseudoscalar and 
vector mesons they constitute which are not tightly con- 
strained by symmetries, the rainbow-ladder truncation of 
QCD's DSEs should produce results that, when measured 
in units of mass, are uniformly « 35% too large. The sys- 
tematic implementation of corrections to this truncation 
then shifts calculated results so that reliable predictions 
and agreement with experiment can subsequently be ob- 
tained. (See, e.g., Refs. 0, [H, Q EE Eft) For each 
value of m, computed observables are approximately con- 
stant on the domain uo = ud(m) ± 5uo, Suo = 0.07; i.e., 
5uo/ud(m) ~ 0.2. 

To define the Faddeev equation we compute the quark 
propagator in Fig. [I] from the following gap equation 

Sip)- 1 = Z 2 (i 1 -p + m bm ) + X(p), 

,A 

£ tu (p) = - / K%(q,p;P)Ssr(q), (4) 

Jq 

where is a Poincare invariant regularisation of the in- 
tegral, with A the regularisation mass-scale, and m bm (A) 
is the current-quark bare mass. Z2(C 2 ,A 2 ) depends ad- 
ditionally on the renormalisation point, £ and the gauge 
parameter. NB. We work in the isospin symmetric limit 
throughout. 

The Faddeev equation requires diquark masses and 
correlation amplitudes. They are obtained from the fol- 
lowing Bethe-Salpeter equation [2[: 

r» P) = J\x qq (i; p)]«- K r s (i, k- p) , (5) 

where: k is the relative and P the total constituent mo- 
mentum; x qq (hP) '■= S(l+)T qq (l;P)S(-l-) T with "T" 
denoting matrix transpose, l± = I ± P/2; and T qq (k;P) 
is the correlation amplitude in a given colour-Dirac chan- 
nel. Equation ([5]) only has solutions in the (3) c channel 
so we can write Y qq = H a T^ P , where the colour is ex- 
pressed through Gell-Mann matrices: {H 1 = zA 7 , H 2 = 
— iA 5 , H 3 = zA 2 }, and T® P is a Dirac- matrix- valued func- 
tion whose explicit form depends on the J p (spin/parity) 
of the correlation. (T qq (l;P) = C^T(-l'P) T C, where 
G = 7274 is the charge conjugation matrix.) 

It has long been known that the lightest diquark cor- 
relations appear in the J p — + , 1 + channels. This can 
be seen to follow from the fact that Y D JPC := Y q * P & 
satisfies exactly the same Bethe-Salpeter equation as the 
J~ p colour-singlet meson but for a halving of the coupling 
strength. Hence only the + and 1 + diquark correlations 
are retained in approximating the quark-quark scattering 
matrix to arrive at the Faddeev equation in Fig.[TJ 

In Fig. [2] we present diquark masses computed with the 
interaction of Eqs. (pQ)-(|3]). The bandwidths show their 
sensitivity to variations in uo. As diquarks do not appear 
in the strong interaction spectrum their masses are not 
observable. We note that the cj-band on m av — m sc is 
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FIG. 2: Thick bands: Evolution with current-quark mass, m, 
of the scalar and axial- vector diquark masses: m sc and m a v 
Bands demarcate sensitivity to the variation uo = Q(rh) zb 5uj: 
lower (upper) edge corresponds to largest (smallest) 00 value. 
(mn, calculated from rainbow- ladder meson Bethe-Salpeter 
equation: rh — 6.1 MeV — 0.138 GeV.) Solid curve: 

Evolution of p- meson mass [111 ]. (NB. The interaction of 
Eqs. JU)-© was deliberately constructed to yield a calcu- 
lated overestimate of m p (m%).) With m p we also depict re- 
sults from simulations of lattice-regularise d Q CD [171 ] along 
with an analysis and chiral extrapolation [18], short dashed 
curve. Thin band: Evolution with m of the nucleon mass 
obtained from our Faddeev equation: rh — 6.1 MeV, Mn — 
1.26(2) GeV. For comparison we provide results from lattice- 
QCD [19|,|20|] and an analysis of such results [21], dashed curve. 



much narrower than that on the individual masses and 
that this difference falls with increasing current-quark 
mass. Since the A-baryon may only involve axial- vector 
diquark correlations, the A-N mass splitting is corre- 
lated with m av — m sc . Hence we infer that Ma — Mn 
will depend weakly on uo and fall with increasing m 2 . 
Notwithstanding the correlation, near agreement be- 
tween the experimental value of Ma — Mn = 0.29 GeV 
and m av — m sc = 0.27(3) GeV at the physical pion mass 
is incidental [22[. (NB. Hereafter parenthesised numbers 
indicate response to variation uo = ud(m) ± 5uo.) 

To complete the Faddeev equation's kernel one must 
specify how the diquark correlations propagate. This in- 
volves a definition of the diquark amplitudes off-shell and 
calculation of the diquark propagators. For the former we 
observe that a diquark's correlation amplitude possesses 
one dominant piece; viz., that single Dirac-amplitude 
which would represent a point particle with the given 
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quantum numbers in a local Lagrangian density, and a 
predetermined number of subdominant amplitudes. In 
taking a diquark off-shell we require that for a correlation 
with large spacelike total momentum; i.e., P 2 ^> M 2 , 
where M is the relevant diquark mass, only the dominant 
amplitude remains pertinent. This is effected by noting 
that in any relevant tensor basis some of the subdomi- 
nant amplitudes involve a factor of P, where P 2 = — 1 
and P is the total momentum of the correlation. We 
multiply each P by a factor h(x) = (l/i)(x/(x + 2)) 1 / 2 , 
where x = P 2 /M 2 . Each subdominant amplitude is sub- 
sequently multiplied overall by a factor g(x) — l/(x + 2), 
while the Lorentz scalar functions associated with each 
matrix-valued tensor are fixed at the on-shell forms ob- 
tained in solving Eq. ©. (NB. On shell, P 2 = —M 2 and 
h(—l) = 1, g{— 1) = 1.) This is a simplification of the 
prescription detailed in Appendix B of Ref. [8[ . 

The propagators for the diquark correlations in Fig.[T] 
can be written in the form (x = P 2 /M 2 ) 



D-\P 2 ) 
(D-V(P 2 ) 



Ml [A sc + p sc f(x) + Q(x)], (6) 
M 2 av [A av S„ v + p av r„ v {x) + Q„ v {x)] . 

(7) 



On-shell the canonical Bethe-Salpeter normalisation con- 
dition ensures that in an internally consistent trunca- 
tion: A = -S (T )(-1), 13 = 1- Q' (T) (-1), .F(-l) = 0, 
^ 7 ( T )(— 1) = 1, where the subscript "(T)" is only relevant 
to the axial-vector propagator and identifies the trans- 
verse piece. The second diagram is denoted by "Q" in 
Eqs. (|6|), ©. It is an easily evaluated, convergent one- 
loop integral whose integrand involves the dressed-quark 
propagator and the appropriate diquark correlation am- 
plitude. The first diagram is associated with T . It rep- 
resents a two-loop integral that must be computed on 
a large domain in the complex-P 2 plane. In principle 
it can be evaluated using Monte-Carlo methods. How- 
ever, we circumvent that additional computer- time con- 
suming step by employing the following par ametrisat ions: 
T(x) = [l+x/(x+2) 3 ]/4; T^{x) = 5^ [l-l/(x+2) 2 ]/2. 
They are: based on an analytic analysis of the integrals; 
correct on-shell; and in qualitative accord with Monte- 
Carlo evaluations at real-P 2 , with fair quantitative agree- 
ment at large spacelike momenta. 

All elements in the Faddeev equation's kernel are now 
defined and computed. (Its explicit form can be con- 
structed following Sec. 5 in Ref. [§] or App. A in Ref. [Io[.) 
The equation is solved using standard methods. Our 
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FIG. 3: Sum of diagrams that determine the diquark prop- 
agators in Fig.[T] T represents the diquark correlation am- 
plitudes; the single line, dressed-quark propagator; and K _1 , 
inverse of the ladder- kernel in Eq. (PQ)-©. 



result is depicted in Fig.O Notably, despite the large 
cj-dependence of the unobservable diquark masses, the 
nucleon mass is only weakly sensitive to this model pa- 
rameter. Our calculation directly correlates an efficacious 
one-parameter model of pseudoscalar and vector mesons 
[7[ with the prediction of a nucleon observable. In par- 
ticular, we have a systematically improvable continuum 
prediction for the evolution of the nucleon mass with a 
quantity that can methodically be connected with the 
current-quark mass in QCD. 

We employed a kernel constructed carefully to ensure 
that corrections to the rainbow-ladder truncatio n sy stem- 
atically move results into line with experiment [11]. The 
corrections can be divided into two classes: nonresonant 
diagrams and meson loops. Both contribute materially at 
small m 2 but vanish with increasing current-quark mass. 
At the physical pion mass we predict Mn = 1.26(2) GeV 
and hence that corrections to the rainbow-ladder trun- 
cation reduce the nucleon mass by 320 MeV. A relevant 
comparison is provided by the analysis of lattice- QCD re- 
sults in Ref. which yields M* =: a N = 1.27(2) GeV. 
Furthermore, it is important and internally consistent 
that at the physical pion mass Mn /M^ pt = 1.34 and 
m p/ m p Xpt = 1-28. Likewise, it is significant that the 
difference between lattice results and our calculations di- 
minishes with increasing current-quark mass. 

Our curve for the nucleon mass appears to lie below 
the lattice results at the largest current-quark masses. 
This might indicate that our Faddeev kernel provides too 
much attraction. That is conceivable given our replace- 
ment of the quark-quark scattering matrix by only the 
lightest correlations in the scalar and axial- vector diquark 
channels. This is a truncation in addition to rainbow- 
ladder and its effect can be quantified. In this connection, 
however, we note that the lattice length-scale (Sommer) 
parameter has decreased with the accessibility of lighter 
current-quarks [24|. Taking that into account for the p- 
meson trajectory would raise all of our calculated curves. 
Therefore in our view Eqs. ([I])-© are satisfactory until 
the Sommer parameter stabilises. 

With the Faddeev amplitude in hand one can calculate 
nucleon form factors as a function of squared momentum 
transfer, Q 2 . The construction of a nucleon-photon ver- 
tex that fulfills the Ward-Takahashi identity for on-shell 
nucleons described by the Faddeev amplitude is detailed 
in Ref. [25[. In this approach one resolves the diquarks' 
substructure. 

Diquark correlations are not point like. Hence, with 
increasing Q 2 , diagrams in which the photon resolves a 
diquark's substructure must be suppressed with respect 
to contributions from diagrams that describe a photon 
interacting with a bystander or exchanged quark. These 
latter are the only hard interactions with dressed-quarks 
allowed in a nucleon. This suppression was not expressed 
in the current of Ref. [8] and that leads to spurious re- 
sults; e.g., the proton electric form factor does not ap- 
proach zero with increasing Q 2 . Here we introduce a 
second, major improvement of Ref. [8]; viz., we refine the 
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nucleons's electromagnetic current and alter as described 
below the expressions in App. D of Re f. 

A key modification is to follow Ref. [11] and write the 
dressed-quark-photon vertex as (y — Q 2 /m 2 p ) 



y 



m p y 



Y T% r^ c (/c; Q) e 



■9(y) 



r»(k;Q)=T» c (k;Q) 

(8) 

Here the first term is the regular Ball-Chiu Ansatz 
[H, [27} . The second introduces a resonant p- meson con- 
tribution, which is explicitly excluded from the impulse 
current when one employs a Ball-Chiu vertex calculated 
from the solution of the rainbow gap equation Q . In this 
piece: f p is the calculated on-shell p-meson leptonic decay 
constant [28[; is a projector transverse to the four- 
vector Q; and Y v vc (k\ Q) is the calculated and canonically 
normalised rainbow-ladder p- meson Bethe-Salpeter am- 
plitude, whose extension off-shell follows the pattern out- 
lined above for diquark correlations. The factor involving 
9(y) — y 2 ) (1 + 2/), with pi 5 2 parameters, accounts 

for any additional effect owing to an off-shell p-meson. 
Emulating Ref. [ll| , p\ was determined by requiring that 
Eq. (|5J) reproduce the DSE 7r-meson charge radius curve 
in Fig. 6 of Ref. [Hj]: Pi = 0.001 + m^/(3.72 GeV) 2 . p 2 is 
relevant to the Q 2 -evolution of form factors. 

To ensure that a diquark's contribution to the nucle- 
ons' form factors is suppressed with increasing Q 2 it is 
necessary to modify the so-called seagull terms in the 
current. If the seagull vertices employed in Ref. [8[ are 
denoted generically by M^ T , then a minimal correction 
is provided by 



M u = M, 



WT 



WT 



(9) 



where m(y) vanishes at y = and as y — > 00 but is oth- 
erwise constrained little. Since a p- meson contribution 
to the quark-photon vertex will be communicated to the 
seagull terms in any internally consistent calculation of 
these five-point Schwinger functions, we employ 



y 



-ps(l+y) 



yj2m p 1 + y 



(10) 



with ps a parameter that governs the degree of suppres- 
sion with increasing Q 2 . m(y) in Eq. ([9]) has no effect on 
the nucleon's static electromagnetic properties. 

Calculation of the nucleons' electromagnetic form fac- 
tors is now a straightforward numerical exercise. The 



TABLE I: Static properties calculated with rh — 6.1 MeV: 
magnetic moments in magnetons denned by calculated nu- 
cleon mass, Mn- Row 2: Results determined by Ref. [Io| ]. 
Row 3: Experimental [30] or inferred [3l| values. 



Calc. 
Ref. [10] 
Expt. 



Pn 



M n r n M M p 



-1.58(3) 2.56(5) 0.00(23) 3.76(5) 3.43(5) 3.47(5) 
-1.51 2.55 -0.60 3.41 3.05 2.99 
-1.91 2.79 -2.63 4.16 4.19 4.07 
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FIG. 4: Shaded band - Calculated result for 
li p G p E (Q 2 ) /G P M (Q 2 ). Bandwidth delimits the sensitivity 
to uj on the domain u) — od(rh) =b 5u: larger (smaller) a;, lower 
(upper) border. Data: open circles, Ref. [32]; filled squares, 
Ref. [33]; filled circles, Ref. [3]; and filled diamonds, Ref. [35| . 



results correlate a one-parameter model of pseudoscalar 
and vector meson properties with the prediction of a 
range of nucleon observables. A raft of results will be 
reported elsewhere. Herein we describe only the static 
electromagnetic properties, listed in Table H and the pro- 
ton form factor ratio, depicted in Fig.[H 

Regarding static properties, the p-meson piece of 
Eq. ([8|) provides ~ 50% of r\, as was also the case 
for the pion [11]. Our kernel deliberately omits meson 
cloud contributions, whose effect was discussed most re- 
cently in Ref. [10]. They contribute uniformly to improv- 
ing agreement with experiment in Table H In accord 
with the model of Ref. [10], we predict quark-core values: 
Pp = 2.36(5), Pp = 0.20(1), where pf means the contri- 
bution from a quark of flavour / to the proton's magnetic 
moment. Charge symmetry entails p^ — —0.40(1), p^ = 
— 1.18(3). NB. At each current-quark mass considered, 
Fig.[2j no calculated magnetic moment or charge radius 
varies by more than 4% on the domain uo = Cd(rh) ± duo. 

Of interest is the behaviour of the neutron's Dirac 
form factor. One can write FJ l (Q 2 ) = (2/3)[/f n (Q 2 ) - 
fi d (Q 2 )}. Our framework predicts F?(Q 2 ) < with an 
evolution at small-Q 2 that can be characterised by two 
Dirac radii; namely, rf n = 0.85(2) fm, r^ d = 0.76(1) fm. 
That our calculated ii-quark contribution to F™ evolves 
more rapidly than that of the d-quark is consistent with 
contempo rary parametrisations of experimental data; 
e.g., Ref. [36]. It owes to the presence of axial- vector 
diquark correlations in our Faddeev amplitude. 

The neutron's Faddeev amplitude involves the quark- 
diquark flavour structures d[ud] + and {^2u(dd)i+ — 
d^d)^}, with the axial- vector correlation being more 
massive than the scalar. The Faddeev equation ensures 
that the two (i-quarks are each equally likely to be found 
in a diquark correlation. Hence, the ud correlations do 
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not favour localisation of the ix-quark at the neutron's 
centre-of-mass any more than they do the <i-quark. How- 
ever, a localisation of the d-quarks does follow from the 
significant probability that the neutron's two d-quarks 
are embedded in a massive 1 + diquark correlation. (NB. 
We have assumed isospin symmetry and hence properties 
of Ff can readily be inferred from this discussion.) 

The proton form factor ratio is depicted in Fig.JH 
Our truncation omits pseudoscalar meson cloud contribu- 
tions. We therefore have r p E > r p M and disagree with ex- 
periment for Q 2 < 2 GeV 2 . (This also explains our inad- 
equate description of r E .) Since pseudoscalar mesons are 
not point like, such contributions rapidly become unim- 
portant beyond this scale and hence the quark core de- 
scribed by our Faddeev equation should be quantita- 
tively reliable. (See Sec. 4.2.3 of Ref. [22[.) We report 
fi p G p E (Q 2 )/G p M (Q 2 ) because its Q 2 > 2 GeV 2 behaviour 
is sensitive to the nucleon-photon current and therefore 
Eqs. (|8l)- (fT0|) become important on that domain. We 
find (^2,^3) = (0.001,0.075) optimises ag reement with 
the polarisation transfer data |33l l34l 1351 ]. Variations 
of ±20% have a noticeable impact but no physically ac- 
ceptable choice for p2,3 can reproduce the Rosenbluth 
data [32]. This is the only parameter fitting herein. NB. 
Calculating the dressed-quark-photon vertex would ob- 
viate the need for the Ansatze and parameters pi,2,3 in 
Eqs. ([8])- ([T0]) and is a natural next step toward true pre- 
dictions of the large-Q 2 behaviour of nucleon form fac- 
tors. 

Our description of diquark propagation is inadequate 
for Q 2 > 4 GeV 2 whereupon its failure to implement di- 
quark confinement leads to pinch singularities within the 
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